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Abstract 

The acoustic spacetime corresponding to perturbed Friedman-Lemaitre-Robert- 
son-Walker universe inherit the space isometries from the original FLRW model, but 
essentially differs in dynamics. The scale factor manifestly depends on the equation of 
state of the matter content. Despite the higher complexity of the background evolution 
the perturbation equation in this space is substantially simpler: the density pertur- 
bations obey d'Alembert equation. Canonical formalism reconstructed in the acoustic 
spacetime enables one to employ the Klein-Gordon scalar product. Consequently, the 
Fourier decomposition of the perturbation field provide the time-independent Fourier 
coefficients and the time-independent spectrum. The perturbation spectrum does not 
depend of the choice of the Cauchy hypersurface from which the data are collected. 
Noether constants associated with the six-parameter isometry group define the com- 
ponents of the momentum, hyperbolic momentum and angular momentum of sound. 

1 Introduction 

Cosmological acoustic field originates due to small perturbation of the Robertson-Walker 
background and contributes to the microwave background temperature fluctuations. The 
perturbation formalisms describing this phenomenon, initiated by Lifshitz [lj, Lifshitz and 
Khalatnikov [2], Landau and Lifshitz [3] and continued by many different authors and groups, 
can be divided into four categories: 

1. Formalisms which adopt gauge- invariant combinations of metric potentials, metric com- 
ponents and their first derivatives. The propagation equations for them are easy to derive 
Bardeen [4], Brandenberger, Kahn and Press [5], but the interpretation of these vari- 
ables is not straightforward Bardeen [4], Katz et al. [6], Deruelle et al. [7j, and Deruelle 
and Uzan [8]. 
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2. Gauge-specific formalisms Q Particularly these, which appeal to the orthogonal gauge jU 
E], the zero shear gauge [Jj, or the longitudinal gauge [TU] . 

3. The gauge invariant formalisms operating with the hydrodynamic quantities like shear, 
divergence of the fluid flow, density gradients or Laplacian [HI [121 H31 CHI EE]. These 
quantities are local, and as such, obviously are gauge- invariant. Their geometrical inter- 
pretation is clear but it is difficult to construct the propagation equations in the partial 
differential form. 

4. The formalisms operating with gauge invariant variables of the desired dynamical prop- 
erties [16j [T7J HH HH]. The theories have Hamiltonian form. Perturbations are identified 
with the scalar fields of time- dependent mass. No particular geometrical interpretation 
is developed. 

More about the mutual relationships between different formalisms can be found in Brechet 
et al. [20] • Each category gives insight into different aspect of the perturbation problem. 
Non of them is known as complete and commonly accepted theory of cosmological density 
perturbations. 

The theory of sound in non static media has developed in roughly the same period. Re- 
search refers both to theoretical [2D [221 1231 122 [221 [2S1 121] and engineering [231 [2S1 ESI ED E2] 
aspects. The key feature that distinguishes the acoustics of fluids in motion is the absence 
of the conserved quantities: energy, momentum and angular momentum [33J of sound. Non- 
conservation of the acoustic energy or momentum leads to the idea of the energy-momentum 
interchange between the wave and the background flow. The more radical formulation states 
that no momentum at all can be assigned to sound wave propagating in non-static environ- 
ment (see Mclntyre |33j critical review). 

Acoustic geometry [3HEH1ES] provides the good insight into the acoustic energy- momentum 
problem and points out the essential difficulties [37] . For arbitrary fluid flow the Unruh met- 
ric [Ml 133 ES] may not admit isometries, and consequently, the energy or momentum cannot 
be introduced as the Noether constants. 

The acoustics of non static media has several tangent points with the cosmological per- 
turbations theory prior the recombination epoch. In both cases: 

a) the environment evolves, and the energy is not conserved, 

b) the density perturbations propagate as waves, 

c) inhomogeneities are small, the linear theory is adequate, 

d) splitting the solutions into the background component and the perturbation is not unique 
(gauge problem). 

Despite these parallels the theory of sound propagation in fluid flows remains unknown 
for cosmologists. To our best knowledge the questions what is the momentum of the sound 
wave [3E] has never been asked in cosmology. 

1 Cosmological literature does not precisely distinguish between the gauge-invariant and gauge-specific 
approaches. There is a tendency to give the name of gauge-invariant to all these theories that eliminate the 
pure gauge perturbations, even if they do eliminate the gauge freedom itself. The expressions such as "gauge 
invariant variables like the curvature perturbation on a constant Hubble surface" are often met and rarely 
criticised. 
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An attempt to answer this question is the aim of this paper. The problem is simpler 
than in classical acoustics, because the Robertson-Walker spacetime admits the six param- 
eter group of isometries [39]. Killing vectors are spacelike. The resulting Noether constants 
can be easily identified with the momentum and angular momentum components. Follow- 
ing Unruh [3H [33] , Visser [36], Golda and Woszczyna [lOl HTJ 02] we construct the acoustic 
spacetime in which the perturbations behave as the scalar field (the propagation equations 
have the d'Alembert equation form). Then, we introduce the energy momentum tensor for 
such field and, following classical methods [SI El El S3] , we construct the conserved currents 
and constants of motion. 

We concentrate on adiabatic perturbations. We do not specify any particular equation of 
state, we claim however that p = p(e) > and 5p/5e > (the speed of sound c s is positive). 
We suppose that some difficulties to express the behaviour of multi component fluids in the 
framework of general relativity (with the famous "frozen amplitudes" problem included) may 
come from the insufficiently precise formulation of the expanding fluid acoustics. The single 
fluid case is good for study the issues. 

We adopt the synchronous system of reference, where the gauge freedom is limited by the 
constraint Sg^^ = 0. This constraint is a good compromise between full gauge freedom and 
completely fixed gauge theory. Computations in synchronous systems are substantially easier, 
while the residual gauge freedom still admits changes of the constant time hypersurfaces. The 
gauge-invariance evoked below means invariance against this residual freedom. The freedom 
is enough to introduce the hyper surf ace-independent variables and the time-independent 
perturbation spectrum^] 

This paper consists of 8 sections. The short section [2] revokes the definition of scalar 
perturbations on the Robertson- Walker background. In the section [3] we derive the partial- 
differential analogue of the Lifshitz-Khalatnikov equations. The gauge-invariant variables are 
introduced in section HI Propagation equation for the new variables is reduced to the wave 
equation form (section [S]). Construction of the acoustic spacetime for the cosmological density 
perturbations is given in section We reduce the propagation equation to the d'Alembert 
equation in this space. Isometries, the momentum of acoustic field and its spectrum are 
discussed the section [71 Some concluding remarks are given in section [HI A short code written 
in Mathematica is attached in the Appendix (section [9]) . 

Throughout this paper we use the convention 8irG = 1, c = 1. The following notation is 
adopted: 

a (v) — conformal scale factor, 
rj — conformal time, 

x k = {x, y, z} — Cartesian 3-space coordinates, 

x M = {rj, x, y, z} — Cartesian spacetime coordinates, 

^9mn — the metric of the maximally symmetric 3-dimensional space. The curvature K is 
an arbitrary real number. 

Where the confusion is unlikely, we abbreviate the notation by writing x for x M and x 
for x k . Dot in x-x stands for x 2 + y 2 + z 2 . We never explicitly express the time or space 

2 We are sure that the Sgo^ — constraint can be relaxed with no harm to the presented results. 
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dependence of the metric tensor g (we hope, this is obvious), but we always keep this depen- 
dence in the metric corrections C(x) and £(x). Vertical lines as in £(x)J J mean the covariant 

differentiation in the space ^g mn . The same £(x)J 3 ' is later abbreviated to A£(x). 

2 Scalar perturbations 

Consider the Robert son- Walker spacetime as the ground state, 

9t*» = a 2 (v) ( q (3)! J > (!) 



mn. 



where ^'g mn stands for the maximally symmetric 3-dimensional space metric tensor, ex- 
pressed in Cartesian coordinates 

(3 W«= (l + ~Kxx^ 5 mn . (2) 

We introduce the small correction Sg^ u 

9w — > 9fMu + Sg^ ( 3 ) 

demanding that 

Sg mn = V m V n £(x) + ~ (C(x) - A £(x)) g mn , (4) 

&9vo = 0. (5) 

V m and A denote the covariant derivative and the Laplacian, respectively — both are calcu- 
lated in the static space % n . C(x) and £(x) are arbitrary functions dependent on space and 
time. The correction (T4]) and (j^J) is the most general formula for the scalar perturbations jH] 
in the synchronous system of reference. 

3 Propagation equations for metric potentials 

The Einstein equations 

Gf_ iu = (6) 
with the hydrodynamic energy momentum tensor 

Tt = diag(-e,p(e),p(e),Ke)) (7) 

are expanded to the first perturbation order with respect to the metrics corrections. The 
zero-order equations reduce to the Friedman equations 

a' 2 (ri) K . . 

3 ^M + 3 ^) = (8) 

a 3 (7?) + a%) a 2 (7?) PW ' 1 J 
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which define the dynamics of the background spacetime. In the first perturbation order one 
obtains the Lifshitz-Khalatnikov equations 



**% + ~ h>\ -hJ j + ^ + 2a - AKhf = 0, for i + j, 



d 2 h , 9n dh 1 . 9w , 
— + (2 + 3c 2 )atf- + -(l + 3c s 2 )(/ i 



m\l 



h/ - 2Kh) = 0. 



(10) 
(11) 



where h is the abbreviation for the metric perturbation: h mn = 5g mn , h = g mn h mn , the 
stroke stands for the covariant space-derivative V n (Tj n = V n T) and c s stands for the sound 
velocity: c 2 = p' '(rj) / V '(77) . i7 = a' {rf) / a 2 {rf) is the Hubble parameter. In this notation the 
scalar perturbations dl])-© read 



hi = -(C(x)-£(xV<)<V+£(x)|^, 
h = hi l = C(x) 



(12) 
(13) 



Inserting f[T2"l) and f[T3"j) into the equations ffTUl) and ffTTl) one obtains the propagation equations 
for the metric potentials C and £ in their partial differential form 



d 2 
drj- 



d 



■£\m n + 2a H—S\ m ^ n + 5 ,m fl ,pi — -£\iplm 9^9 

1 



ip An 



+ £{*)\ imp i{9 ll 9 np - 9 in 9 pl ) ~ ^C\J n - 4KS\J n = 0, for m + n, 

— C + (2 + 3c 2 )atf-C + (l + 3c 2 ) 



x 



\s\i hsm {g lm g hs - \g lh g sm ) - \ {C/ + 3KC) 



o. 



By adopting the identity 

£\ ipm i(g m g pl - g rp g nl ) + ^ mpl {g u g np - g m g pl ) - 4K£ { J n = o, 

one can reduce the first equation f[T4"|) to 



^8 + 2aH^-£+ 1 -(8/-C) 
or} 1 or/ 3 



n 



0, for m ytz n. 



- \m 



(14) 



(15) 



(16) 



(17) 



To fulfil the equality f(x^)\ rn ) n = 0, Vm ^ n in arbitrary coordinate system, the function 
f{x^) must be a constant. Then, / obviously satisfies f(x^)\rn Tn — which implies also that 



r5 2 (9 1, 

wl S + 2aH— £ + -(£/- C) 
or\ A or/ 3 



0. 
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On strength of the identity 

C f,,lm„sh Ah „sm\ OZSC \ m n 

0\lhsm{9 9 ~9 9 ) - 2A<t| m ' = (J, 

the second equation ([15]) simplifies to 

+ (2 + 3cl)aH^-C + ~(1 + 3c s 2 ) [[£,,1' - C] ,J m + 3^/ - C 

Finally, the propagation equations take the form 
d 2 AS a' dAS 1 



dr] 2 
d 2 C 



a dr] 



A (AS - C) , 



dr] 2 



-(3c s 2 + 2) 



a dr) 



+ c 2 s (A + 3K) (AS - C) 



(19) 
(20) 

(21) 
(22) 



These equations are the partial differential analogue to the Lifshitz-Khalatnikov system. 
Neither the metric potentials nor the density contrast evaluated from them 



fe(x) 



3a 2 



Z-— + (A + m)(AS-C 
a or] 



(23) 



are observables. All of them are ambiguous due to the existing gauge freedom. 



4 Gauge-independent variables 

When the gauge is fixed, like the fluid orthogonal gauge or the zero shear gauge, all the 
quantities are gauge-specific. There is no standard technique to introduce gauge-invariant 
quantities on fixed hypersurfaces, neither to compare gauge-specific quantities between dif- 
ferent gauges. For the properly constructed gauge-invariant spectrum the existence of gauge 
freedom is indispensable. 

In the case we discuss in this paper, the gauge freedom is guaranteed by two arbitrary 
space dependent functions Qi(x k ) and £2(2^) • The gauge solutions for AS and C are 



AS g (x») = -g 1 (x k )-2 ^d v AQ 2 (x k ), (24) 

J a (v) 

c g (x») = gi ( x k ) + 6 ^p\g 2 (x k ) + 2 f-Ldq Ag 2 (x k ). (25) 

a 2 (r]) J a(rj) 

Not eliminating this freedom, we look for quantities of physical interest, which are indepen- 
dent of Q\(x k ) and Q 2 (x k ). 

Below we abbreviate the notation by writing x for x^ and x for x k . Our construction 
of the gauge-independent variables is based on the fact that the pure-gauge energy density 
perturbation 

5e 9 (x) = &(x)^ = &(x)e(t) (26) 
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forms the product of 02 ( x ) an d the time derivative of the background energy density e(t) 
(the formula dr] = a(t)dt connects the two different time parameters). The formula (126 p can 
be easily obtained by inserting the system f )2^|) and (125]) into ( 123]) . Fictitious perturbations 
for other scalar quantities (the expansion rate, the scalar curvature, etc. ) have the same 
structure 





= 02 M, 


"50(x)~ 


= 02 (X), 


"5i2(x)" 


l m . 
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02 (x), etc. 



(27) 



Differentiation with respect to time removes the gauge dependence Qi (x). Hence, a simple 
mnemotechnic rule arises, that allows one to convert any fractional perturbation of scalar 
quantity (the so called contrast) into the corresponding dimensionless gauge-independent 
variable 



MX)" 



59 (x) 



'59 (x) 



L J ' 



5i?(x) 



'5i2(x) 



Li?(t) 



etc. 



(28) 



The rule consists in completing the expression for contrast with two derivative dots: one in 
the denominator, and one over the whole fraction. 



5 Propagation of sound 

Arbitrary time dependent factor f(t) multiplying the perturbation amplitude K(x) 

r(x) — ► /(t)y(x) (29) 

does not affect the gauge-invariance. The formulae ( 12 8 p define the families of gauge inde- 
pendent quantities, with the amplitude scaling freedom ( l2"9~j) . The choice of the factor f(t) is 
one of the possible tools to find the canonical form for the propagation equation. We guess 
the factor / in the form f(t) = a 2 (t), and introduce 



7(x) 



a\t) 



(30) 



to define the gauge-invariant variable for the energy density perturbation. Returning to the 
conformal time 77 we write 



<9 2 7 (x) 
dr\ 2 



a(rj) c s (r])\ dr] 



c s 2 (r / )A 7 (x) = 0, 



(31) 



where the function a(rj) reads 



0(77) = a(rj), 



pjrj) + ejjj) 



(32) 
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The background dependent quantities: Hubble parameter H, pressure p, the energy density 
e and the sound velocity c s involved in ( ]52|) are uniquely determined by Friedman equations 
and the equation of state. 

The equation ( 15Tj) describes the acoustic waves propagating in the expanding environment 
with variable sound velocity c s (rj). 

Proof of the equation f |5T|) consists in simultaneous computer aided reduction of the 
system of equations: <gl$, $2$, A snort code written in Mathematica 

is presented in the Appendix. We use normal forms of commands and expressions, where 
possible. The full code for sections [2H5] written in Tensorial package is available from [4*5"] . 



6 Acoustic spacetime 

Appropriate redefinition of the time parameter completes the construction of the acoustic 
spacetime [MJ [35j ES]- To avoid the explicit time- dependence of the sound velocity c s (ry) in 
the equation ( 15Tj) we introduce 

C= fcMdv. (33) 



The resulting propagation equation 

<9 2 7(?) , a'(C)97(?) 



A 7 (y) = 0, (34) 



d( 2 o(C) d( 
is the explicit form of the d'Alembert equation 

fiWVVTfr) = 0, (35) 
in the spacetime j: = {(, x, y, z} with the metric 



° 2 (0 ( n 1 (3)1! ) > where a (v(0) = a \ W^lH- (36) 



^ - w ^ v (3 W ' v " V 3^ 

Equations (EIJ) and (155]) say that the sound in the Robertson-Walker spacetime M. with 
the scale factor a(rj), propagate as massless scalar field (1551) in another Robertson- Walker 
spacetime M.' (1361) . Sachs and Wolfe [16] reported wave-like behaviour of scalar perturbations 
in the flat, radiation-dominated universe. The equations (135|) and (1361) may be considered as 
the generalization of the Sachs- Wolfe theorem to spatially curved Robert son- Walker models 
and to arbitrary equation of state of the barotropic form p = p(e). 

In the following sections we drop the sign prime. The coordinates {(, x, y, z} we briefly 
denote by the same letter x. 
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7 Symmetries, momenta and spectrum of sound 



We introduce Lagrangian for the acoustic field 

1 
2 

and derive the canonical energy-momentum tensor 



£ = -7^<9" 7 (x)<9 CT 7(x) (37) 



T» = 0" 7 (x) 0" 7 (x) - -Q^QpcdPjix) «9 CT 7 (x). (38) 

The propagation equation (13"5]) is the Lagrange equation for £, or equivalently, it follows from 
the V V T^ U = identity. Both C and are gauge- invariant, therefore, all the conserved 
quantities constructed by means of them have the same property. The procedure below 
follows that of Katz et al. [6J, Deruelle et al. [7], and Deruelle and Uzan [8] except for 
the background space is (136]) and the field 7(x) is already gauge-invariant. We consider 
conserved currents jf 1 = T llu K^ v . The six space- like isometries immediately provide 
the six Noether integrals tt(K^) 

tc(K®)= [ J®dE' t = 1% V K^ U dE" = const (i) (39) 

J E J E 



where is an arbitrary Cauchy surface. The algebra of Killing vectors, and consequently, 
the interpretation of constants f l39]) depend on the sign of curvature. For Killing basis chosen 

as 



= {K$,K®} 



3 



3 



x-x5 ij -2x i x j - VkYI 



€ ijk X k 



fc=l 



(40) 



Kf j = ^e ljk x k (41) 
fc=i 

K (i)o = K ®o = Q (42) 

generators of infinitesimal isometries 

T^ = -iK^ j d jt L®=-iK$ )j d, (43) 
satisfy the commutation relations 

3 

[ T « ;T (i)] = iv^K^ e ^T (fc) (44) 

k=l 

3 

[L®,L&] = -iJ2e ijk L {k) (45) 
fc=i 

3 

[L (4) ,T W ] = -i^ e y fc T (fc) (46) 



fc=i 
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For K > 0, both and assure the angular momenta conservation (six constants 
of motion). For K = the right-hand side of (j44"l) vanishes and operators T^> generate 
translations. Translations Kj) conserve the momentum (3 constants of motion), while the 
rotations K^ 1 conserve the angular momentum (the next 3 constants of motion). For K < 
vectors Kj) correspond to hyperbolic momentum [l6j STJ HH]. The absence of the time 
isometrjjf] breaks down the energy conservation. 
The Casimir operator 

3 

^T (i) T (i) = -A (47) 
i=i 

is equal to minus Laplacian. T"W and A commutate 

[T«,A]=0 (48) 

hence, each pair {T«,A}, i = 1,...,3 has common eigenf unctions. There are no common 
eigenfunctions for pairs {T® } TV>} with i ^ j and K ^ 0. Solutions u k {x^ ,n l ) to equation 
(I35p which lie in the kernel of operator niU- l \ and are eigenfunctions of the operator n^TW 

mL^Uk^irt) =0 (49) 
ntTWukix^n*) = {k - W^K)u k (xP, rt) (50) 

define plane waves of the wavenumber k, propagating in the direction ri i . On the strength 
of (l4*8|) the same functions are also the eigenfunctions of the Laplace operator 

A« fc (x", rj) = ~(k 2 - K)u k (x", rj) (51) 

Since A is Hermitian, the functions Uk{x^,n l ) of different k are orthogonal on constant time 
hypersurfaces. When separated they read 

u k = ^- X k(OF(x,n,k). (52) 

a (C) 

The time dependent factor Xk(C) obey the evolution equation 

- ^ + ^crao- (53) 

and implicitly depend on the equation of state. The space-dependent solutions F(x, n, k) to 
the Helmholtz equation ()5ip are 



1 f-y/=Kn-x + 1 - fx-x 
(27r")V2 i 1 + fx-x 



H*, n, k) = — t— V J 4 ^ • (54) 



3 Conformal Killing vectors in both spaces Ai and M.' are not identical. The conserved quantities related 
to conformal isometries [51 [71 [S] of the space A4' are a separate issue, not discussed in this paper. 
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In the K = limit F(x, n, k) tend to e lfc ("' x ). For that case the spectral analysis of the acoustic 
field resembles that for electromagnetic field |1HJ EDI EI]- For the negative space curvature 
functions F(x, n, k) are Shapiro functions [52]. Principal series characterised by positive 
wavenumbers k > consists of functions orthogonal and complete in L? [531 El]- Supple- 
mentary series of regular, bounded, non-oscillating and non-orthogonal functions -F(x, n, k) 
with imaginary wavenumber k G (0, \/K) are redundant to expand the square integrable 
perturbations. However, they contribute to Fourier decomposition of weakly homogeneous 
stochastic processes [22 EE]- A single function of supplementary series is also needed for 
completeness in the quantum theory of charge [57]. In both cases, flat (K = 0) or open 
(K < 0) universes, the Fourier bases are denumerable, and the spectrum is continuous. 

For positive curvature the functions F(x, n, k) are known as Sherman- Volobuyev func- 
tions [5H1 EH]- The spectrum of Beltrami-Laplace operator is numerable. The Sherman- 
Volobuyev functions do not form the complete set, neither are orthogonal [BO], and therefore, 
the case of closed universe requires a separate treatment. 

In what follows, we limit ourselves to non-positive curvature and to positive wave num- 
bers (continuos spectrum and principal series). We consider the acoustic wave 7(x M ,n*) 
propagating in the direction n l , i. e. an arbitrary solution to the equation f[3"5"]) that can be 
expanded 



n 



(a k u k (x fl , n l ) + aluKx^, n l ))dk 



in orthogonal basis u k (x^ 



n 



k 



k-i^K 



5(k - k') 



with Fourier coefficients 



a fc = ( 1 - j:^-^ ) ( u k, tW 



Symbol (,)jcg means the Klein-Gordon scalar product 



h,<P2)icg 



(55) 



(56) 



(57) 



(58) 



where stands for Wronskian. For arbitrary pair <\>\ and 02 of complex solutions to the 
equation ( 135]) the divergence of Wronskian vanishes 



n,<p 2 ) 



^2 







(59) 



This assures that the integral (158]) is independent of the choice of the E hypersurface. In 
particular, the integral (158]) and all other quantities defined on this base are invariant under 
the gauge modifications. Coefficients a k are constant in time. Normalisation (156]) together 
with separation ( 15 2 j) result in 

W ll ( X k,ti) = -i (60) 
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and 



Jf(x, n, k)F*(x, n, k') 




(61) 



To express the momentum (l3T?j) of the wave moving in the direction n l we write 



= [ AE°d ol {x», n i )n i K®"d v y(x*, n*) 



(62) 



On the strength of f )55|) . f )56|) . fl57|) and fl60l) after laborious but straightforward calculations 
one obtains 



The quantity Pk = Q-k a*k is constant in time and invariant under both: the gauge transfor- 
mations of the reference system, and the unitary transformations of Fourier bases Mfc(x M , n l ). 
Pk is an observable. It is the intrinsic property of the acoustic field and, according to (|63|) . 
can be regarded as the momentum spectrum or hyperbolic momentum spectrum of sound. 

8 Summary 

The presented gauge- invariant formalism is the first formalism which essentially unifies both 
subhorizon and superhorizon regimes. There is no need to divide solutions into classes accord- 
ing to their length sea This particular property is due to the fact that the perturbation 
variables Sp/p defined in the original spacetime form the massless scalar field 7 in the cor- 
responding acoustic spacetime. Fourier coefficients are based on the Klein-Gordon scalar 
product (|58|) . and therefore, they are gauge-independent and constant in time. Frequency 
distribution of the momentum is described by the hypersurface- independent spectrum Pk. 

The acoustic spacetime enable one to construct the energy momentum tensor of sound. 
The six-parameter isometry group acting in this space results in the conservation of mo- 
mentum, hyperbolic momentum, and angular momentum of sound waves. The momentum 
conservation fl39|) show that the physical nature of waves, whether "inside" or "outside" the 
particle horizon, is the same. The wave momentum does not change at the horizon crossing, 
consequently, the phase and group velocity also remains the same. The evidence of the mo- 
mentum conservation testify that the expected freezing of waves larger than the horizon-size 
is rather a mathematical effect caused by neglecting some terms in differential equations. 

4 Except for non vanishing curvature or cosmological constant. 



7T 




(63) 
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Although the total momentum density in the Robert son- Walker universe is expected to 
be small, the momenta of the constituent waves may be arbitrary. There are no adequate 
physical conditions necessary for exciting and maintaining standing waves [66]. The generic 
solution to the perturbation equation is the running wave. The absence of time-like isometry 
results in the absence of the energy conservation. Acoustic waves may undergo amplification 
and backscattering, due to the same mechanism which works in the gravitational waves 
theory jSHE2]. 

The perturbation spectrum depends on curvature. In the open universe scalar pertur- 
bations are composed of acoustic modes with real wavenumbers < k 2 (principal series) 
and non- acoustic modes with imaginary wavenumbers — 1 < k 2 < (supplementary series). 
Both series of modes may affect the temperature fluctuations, and particularly contribute 
to the lowest multipoles. The dispersion on curvature, dispersion on cosmological constant 
and the gravitational waves amplification affect the same range of scales. This may impede 
to separate these phenomena in the CMBR data. 

The lack of appropriate identification of additive constants of motion in cosmology was 
the essential obstacle to extend the statistical physics theory to the expanding universe. 
Removing this obstacle by reformulating the problem in terms of acoustic geometry enables 
one to replicate the fluid statistical mechanics [67J to universes with Robertson-Walker 
symmetry. Moreover, for the radiation-dominated epoch the resulting acoustic space time is 
static, and consequently, in this particular but important case also the energy conservation is 
restored. The construction of statistical fluid cosmology is beyond the scope of this paper, but 
we appreciate its importance. The classical fluid fluctuations may admix to the temperature 
spectrum, preventing us from knowing which of the observed fluctuations are credible tracers 
of the early quantum processes. 

9 Mathematica 7.0 code — Appendix 

(* Preliminary definitions *) 

Clear All [" Global'*"] 

X = {r],x, y, z}; x = Sequence@@A; 

X = {r]_, x_, y_, z_}; x = Sequence@@A; 

Y = RestLY]; y = Sequence@@y; 
P = {pl,p2,p3}; 

P = {pi-, p2_, p3_}; p = Sequence@@P; 
spacedif = Sequence@@Thread[{Y, P}]; 

zero = Sequence@@Constant Array [0, Length [Y]]; 

(* Space metric tensor *) 

g = DiagonalMatrix[{g[y],fl[y],fl[y]}]; 
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(* Friedman equations *) 

FR2 = p[rj] == -2^jpr + ^pr - ; 

(* Background dynamics *) 

a'[ V _] = a[rj] 2 H[rj}; 

e'[?7_] = -3a [77] # [77] (e [77] +p[rj\); 

H>[ V _] = -a[ V ] (ffW 2 + J(cW+3pW)); 

eps = Solve[Fm,e[?7]][l,l]; 

e[?7_] = Simplify [e [7/] /.eps]; 



(* generalised Lifshitz equations: A = AS, /i = —C *) 
A (2,zero)y = - 2 ^d v X[x] - i(Lap[A][x] + Lap[/j] [a;]); 



// (2,zero)) [ ^ ] = _ ( 2 + 3c[ V ] 2 ) ^d^x] + § (1 + 3cfo] s 



(Lap[A][x] + Lap[/i][x] + 3X(A[x] - ^[x])); 



(* The density contrast *) 

<%] = (-(Lap[A][x] +LapH[x] +3X(A[x] + 3^9^]); 

(* Derivative reccurence commands *) 
A (kJnteger/;fe>2,zero)y ;= d {rhk „ 2} \^ zeT0 ^x] 
// (kJnteger/; fe >2,zero)y ;= g^^^zero)^ 

A (2 ^k]:=VcedifA (2 ' Zer0) W 
^ 2 ^[-]-^pacedi^ (2 ' Zer0) W 

(* Laplacian *) 

Laplasjan[x_, gJ\ := Module[{}, dim = Lengthfx]; 

dim dim , x 

A dim := ^/Detb] 5 £ 9xW (y De %] Inverse [0] [i, fc]d*[ fc ]#lj &] 

Laplasjan[Y,#]; Lap[/_][x] = A dim [/[a:]] //Simplify; 

(* Canonical variable *) 
Clear [H] 
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j[x] = a[ V ]H[ V ]%^//Sim^fy; 
(* The wave equation *) 

wave = d {rh2}1 [x] + [2^1 - fjl] d vl [x] - c^MlM == 0; 
Timing [Reduce [wave /.eps]] 
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